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Abstract 

We systematically investigate the pion transverse momentum spectrum, 
elliptic flow/, and Hanbury-Brow/n-Tw/iss (HBT) interferometry in the granular 
source model of quark-gluon plasma droplets in ultra-relativistic heavy ion 
collisions. The granular source model can w/ell reproduce the experimental 
results of the Au-Au collisions at = 200 GeV and the Pb-Pb collisions at 

^Sivjv = 2.76 TeV with different centralities. We examine the parameters of 
the granular source models with an uniform and Woods-Saxon initial energy 
distributions in a droplet. The parameters exhibit certain regularities for 
collision centrality and energy. 


1 Introduction 

Single particle transverse momentum spectrum, elliptic flow, and Hanbury-Brown- 
Twiss (HBT) interferometry are three important final particle observables in high 
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energy heavy ion collisions. They reflect the characters of the particle-emitting 
sources in different aspects and at different stages. Therefore, a combined investi¬ 
gation of these observables can provide very strong constrains for source models. 

So far, much progress has been made in understanding the experimental data 
of the heavy ion collisions at the top energies of the Relativistic Heavy Ion Collider 
(RHIC) [TJ|2ll3]Il|5ll6]|Zll[9]. However, more detailed investigations of the 
physics beneath the data through multiobservable analyses are still needed. On the 
other hand, the experimental data of the Pb-Pb collisions at ^/sNN = 2.76 TeV at 
the Large Hadron Collider (LHC) have been recently published )10l ITT] IT^ 113] . It 
is an ambitious goal for models to explain the experimental data of particle spectra, 
elliptic flow, and HBT interferometry in different centrality region consistently for 
the heavy ion collisions at the RHIC and LHC. 

In Refs. |141 [TSl fT6]. W. N. Zhang et al. proposed and developed a granular 
source model of quark-gluon plasma (QGP) droplets to explain the HBT data of 
the RHIC experiments [TT] [T^ [THl US] . In Ref. [211122] . the granular source model 
was used to explain the pion transverse momentum spectrum and HBT data of 
the most central heavy ion collisions at the RHIC and LHC. Motivated by these 
successes, we systematically investigate the pion transverse momentum spectrum, 
elliptic flow, and HBT interferometry for the granular sources in the heavy ion 
collisions at the RHIC and LHC energies with different centralities. The granular 
source parameters for an uniform and Woods-Saxon initial energy distributions in a 
droplet are examined and compared. 


2 Granular Source Model 


The granular source model of QGP droplets regards the whole source evolution as 
the superposition of the evolutions of many QGP droplets. Each droplet has a 
position-dependent initial velocity and evolves hydrodynamically. The model con¬ 
struction is based on the following suggestions. In the heavy ion collisions at top 
RHIG energies and LHG energies, the created strong coupled QGP (sQGP) systems 
at central rapidity region may reach local equilibrium at a very short time, and then 
expand rapidly along the beam direction (z-axis). Because the local equilibrium sys¬ 
tem is not uniform in the transverse plane [x-y plane) [23]. The system may form 
many tubes along the beam direction during the subsequent fast longitudinal ex¬ 
pansion and finally fragment into many QGP droplets with the effects of “sausage” 
instability, surface tension, and bulk viscosity [151 [24l [2^ I26j . 

As in Ref. |15| . we suppose the QGP droplets of the granular source initially 
distribute within a cylinder along z-axis by 


dNd 

dxodyodzo 


x0(n,-\zo\). 


9{JZt — po) 


( 1 ) 


Here po = ^^Xq + y^ and zq are the initial transverse and longitudinal coordinates 
of the droplet centers. The parameters TZt and TZz describe the initial transverse 
and longitudinal sizes of the source, and ATZt is a transverse shell parameter (15] . 
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In Ref. [^, the Bjorken hypothesis |27j of longitudinal boost-invariant is used 
to describe the longitudinal velocity of droplet for the most central collisions, and 
the transverse velocity of droplet has a form of exponential power. Considering 
the longitudinal velocity of droplet varying with collision centrality, we introduce 
here also a longitudinal power parameter, which will be determined by experimental 
data, to describe the longitudinal velocity phenomenologically. The initial velocities 
of the droplets in granular source frame are assumed as m 

Vdi = sign(roi) • a, (, f = 1, 2, 3, (2) 


where rpi is xq, yo, or zq for f = 1, 2, or 3, and sign(roi) denotes the signal of 
roi, which ensures an outward droplet velocity. In Eq. Q, TZi = (TZttTZtjT^z), 
di = icLx,ay,az) and bi = {bx,by,bz) are the magnitude and exponent parameters 
in X, y, and z directions, which are associated with the early thermalization and 
pressure gradients of the system at the breakup time. It is also convenient to use 
the equivalent parameters ax = (ax + ay)(2 and Aqt = ax — ay instead of Ux and 
Qy. The parameters ax and Aax describe the transverse expansion and asymmetric 
dynamical behavior of the system at the breakup time, respectively. For simplicity, 
we take bx = by = bx in calculations. The parameters bx and b^ describe the co¬ 
ordinate dependence of exponential power in transverse and longitudinal directions. 

In the calculations of the hydrodynamical evolution of the droplet, we use the 
equation of state (EOS) of the S95p-PCE165-vO |28| . which combines the lattice 
QCD data at high temperature with the hadron resonance gas at low temperature. 
We assume systems fragment when reaching a certain local energy density cq, and 
take the initial energy density of the droplets to be 2.2 GeV/fm^ for all considered 
collisions for simplicity |22j . The initial droplet radius tq is supposed satisfying a 
Gaussian distribution with the standard deviation tr^ = 2.5 fm in the droplet local 
frame |22j . We consider an uniform initial energy distribution in droplet and a 
Woods-Saxon distribution, 

e(r) = ep-^ -, (3) 

-Fl 


where a = O.lrg. 

With the evolution of the hot droplets, the final pions freeze out at temperature 
Tf with the momenta obeying Bose-Einstein distribution. To include the resonance 
decayed pions emitted later as well as the directly produced pions at chemical freeze 
out earlier, a wide region of is considered with the probability |22| 


dP 


_ ^chem 

oc /dire +(l-/dir) 


_ ^chem 

xe , (Tchem > Tf > 80 MeV), 


( 4 ) 


where /dir is the fraction of the direct emission around the chemical freeze out 
temperature Tchem. ATdir and ATdec are the temperature widths for the direct and 
decay emissions, respectively. In the calculations, we take /dir = 0.75, ATdir = 10 
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MeV, and ATdec = 90 MeV as in Ref. j22j . The value of Tchem is taken to be 165 
MeV as it be taken in the S95p-PCE165-vO EOS (28) . 

After fixing the parameters used in the calculations of hydrodynamical evolution 
and freeze-out temperature, the free model parameters are the three source geome¬ 
try parameters {Rt, ^Rt< Rz) and the five droplet velocity parameters (a-r, 

Oz, bx, bz). They are associated with the initial sizes, expansion, and directional 
asymmetry of system, and have significant influence on the observables of pion 
momentum sprectra, elliptic flow, and HBT radii in the granular source model. In 
next section, we will determine these parameters by the experimental data of these 
observable, and examine their variations with collision energy and centrality for the 
heavy ion collisions at the RHIC and LHC. 


Results of Pion Momentum Spectrum, Ellip¬ 
tic Flow and Interferometry 


> 

CD 

O 




O PHENIX, 0-5% centr 
□ PHENIX, 10-20% centr 
O PHENIX, 30-50% centr 
STAR, 0-5% centr 
STAR, 10-20% centr 
STAR, 30-50% centr 
ALICE, 10-20% centr 
ALICE, 40-50% centr 
I I I I ■ I 


Uniform E(r) 
Woods-Saxon ^(r) ' 


1 1.5 


2.5 


Pj (GeV/c) 


Figure 1: (Color online) The pion transverse momentum spectra of the granular 
sources for the RHIC Au-Au collisions at s/snn = 200 GeV and the LHC Pb-Pb 
collisions at =2.76 TeV, for the uniform and Woods-Saxon initial energy 

distributions in a droplet. The experimental data of PHENIX [3T], STAR [52] . 
and ALICE [TU] are also plotted. 


In high energy heavy ion collisions, the invariant momentum distribution of final 
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particles can be written in the form of a Fourier series |29l 130], 


E 


(PN 

d?p 


1 <fN 

2-k pTdpxdy 


1 + ^ 2vn cos(n^) , 

n 


( 5 ) 


where E is the energy of the particle, px is the transverse momentum, y is the 
rapidity, and (j> is the azimuthal angle with respect to the reaction plane. In Eq. 
(§, the first term on right is the transverse momentum spectrum in the rapidity 
region dy, and the second harmonic coefFicient V 2 in the summation is called elliptic 
flow. 

In Fig. we plot the pion transverse momentum spectra of the granular sources 
with the uniform and Woods-Saxon initial energy distributions in a droplet. The 
experimental data of the Au-Au collisions at s/snn = 200 GeV at the RHIC j31l[32] 
and the Pb-Pb collisions at ^/sNN = 2.76 TeV at the LHC j lO j are also plotted. 
In Fig. 1^ we plot the pion elliptic flow results of the granular sources with the 
uniform and Woods-Saxon initial energy distributions, and the experimental data of 
the Au-Au collisions (33] and the Pb-Pb collisions |llj . 



Figure 2: (Color online) The pion elliptic flow of the granular sources for the 
RHIC Au-Au collisions at = 200 GeV and the LHC Pb-Pb collisions at 

y/spfN = 2.76 TeV, for the uniform and Woods-Saxon initial energy distributions 
in a droplet. The experimental data of STAR [33], and ALICE m are also 
plotted. 


The transverse momentum spectrum and elliptic flow of the granular sources are 
well in agreement with the experimental data, except for the elliptic flow results at 
px > 2 GeV/c. The differences between the results of the granular sources with the 
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uniform and Woods-Saxon initial energy distributions are small. The experimental 
data of the momentum spectrum and elliptic flow at the same centralities can si¬ 
multaneously give the strong constraints to the velocity parameters of the granular 
sources. After then, the geometry parameters of the granular sources for the colli¬ 
sions with the different centralities can be further determined by the experimental 
data of HBT interferometry at the same centralities. 
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Figure 3: (Color online) The HBT results of the granular sources for the uniform 
and Woods-Saxon initial energy distribution, and the experimental data for the 
Au-Au collisions at the RHIC [5D] and the Pb-Pb collisions at the LHC [T^ 
with different centralities. 


Two-particle HBT correlation function is defined as the ratio of the two-particle 
momentum spectrum P(pi,p 2 ) to the product of two single-particle momentum 
spectra P(pi)P(p 2 ). It has been widely used to extract the space-time geometry, 
dynamic and coherence information of the particle-emitting source in high energy 
heavy ion collisions j34l l3^ [36l 1371 138] . In the usual HBT analysis in high energy 
heavy ion collisions, the two-pion correlation functions are fitted by the Gaussian 
parameterized formula 


C'(q'out,9side,®Ong) = l+ Ae 


( 6 ) 


where gout, Qside, and giong are the Bertsch-Pratt variables |39][40], which denote 
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the components of the relative momentum q = Pi — P 2 in transverse out and 
side directions and in longitudinal direction, respectively. In Eq. ([^ A is chaoticity 
parameter of source, Rout, Rside, and i?iong are the HBT radii in out, side, and 
long directions. 

We plot in Fig. [^the two-pion HBT results for the granular sources and the 
experimental data of the RHIC Au-Au collisions |20] and the LHC Pb-Pb collisions 
|13| with the same centralities as the experimental data of the spectrum and elliptic 
flow. Here, kr is the transverse momentum of pion pair. One can see that the 
granular source models can well reproduce the experimental HBT radii and their 
variations with kx- The results of the chaotic parameters A of the granular sources 
are larger than the experimental data, because there are many effects in experiments 
can decrease A [Ml [35] |36] [371 138] , which exceed our considerations. 

In Table 1, we present the parameters of the granular sources with uniform 
initial energy distribution in a droplet, determined by the experimental data of 
the momentum spectra, elliptic flow, and HBT radii. The values of the source 
geometry parameters indicate that the sources are a short cylinder with small shell 
effect (small ARt). For certain collision energy, the source geometry parameters 
Rt, ARt and Rz increase with the collision centralities. However, for the 10-20% 
centrality, these geometry parameters for the LHC collisions are larger than those for 
the RHIC collisions. The large difference between the droplet velocity parameters bx 
and bz indicates the different dynamical behaviors of the sources in the transverse 
and long directions. The parameter Aax increases with decreasing centrality. And, 
the values of ax and az are almost independent of collision centrality. A detail 
analysis on the relationships between the source parameters and the granular source 
space-time evolution can be seen in Ref. m- 


Table 1: The geometry parameters (in fm unit) and velocity parameters of the 
granular sources with uniform initial energy distribution in a droplet. 


Centrality 

TZx 

ATZx 

Tlz 

ax 

Aax 

az 

bT 

bz 

RHIC, 0-5 % 

5.8 

0.7 

3.9 

0.469 

0.066 

0.593 

0.76 

0.13 

RHIC, 10-20% 

4.5 

0.5 

2.9 

0.454 

0.115 

0.593 

0.56 

0.11 

RHIC, 30-50% 

2.5 

0.3 

0.5 

0.437 

0.156 

0.593 

0.37 

0.06 

LHC, 10-20% 

6.0 

0.9 

5.5 

0.431 

0.092 

0.592 

0.35 

0.13 

LHC, 40-50% 

2.5 

0.4 

1.8 

0.407 

0.131 

0.590 

0.23 

0.03 


In Table 2, we present the parameters of the granular sources with Woods-Saxon 
distribution of initial energy in a droplet, determined by the experimental data of 
the momentum spectra, elliptic flow, and HBT radii. By Comparing the two set 
parameters in Table 1 and Table 2, one can see that there are some differences 
between the Rz and Uz values for the RHIC collisions. Also, for the LHC collisions, 
the values of ax for for the granular sources with the Woods-Saxon distribution are 
larger. 
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Table 2: The geometry parameters (in fm unit) and velocity parameters of the 
granular sources with Woods-Saxon initial energy distribution in a droplet. 


Centrality 

'R-x 

AT^j- 


dj' 

Acf.j' 

Gz 

b'j' 

hz 

RHIC, 0-5 % 

5.8 

0.7 

5.1 

0.469 

0.066 

0.52 

0.76 

0.13 

RHIC, 10-20% 

4.5 

0.5 

4.0 

0.457 

0.122 

0.52 

0.56 

0.11 

RHIC, 30-50% 

2.8 

0.3 

1.8 

0.453 

0.156 

0.52 

0.37 

0.06 

LHC, 10-20% 

6.0 

0.9 

5.5 

0.496 

0.092 

0.59 

0.43 

0.13 

LHC, 40-50% 

2.5 

0.4 

1.8 

0.434 

0.127 

0.59 

0.23 

0.03 


4 Summary and Conclusions 

We systemically investigate the pion transverse momentum spectrum, elliptic flow, 
and HBT interferometry in the granular source model for the heavy ion collisions 
at the RHIC highest energy and the LHC energy. The centrality dependence of the 
observables at the two energies are examined. By comparing the granular source 
results with the experimental data of the Au-Au collisions at -y/SNN — 200 GeV 
at the RHIC and the Pb-Pb collisions at ^/snn = 2.76 TeV at the LHC with 
different collision centralities, we investigate the geometry and velocity parameters 
in the granular source models with an uniform and Woods-Saxon initial energy 
distributions in a droplet. The parameters as a function of collision centrality and 
energy are examined. Our investigations indicate that the granular source model 
can well reproduce the experimental data of pion transverse momentum spectra, 
elliptic flow, and HBT radii in the Au-Au collisions at = 200 GeV with 0- 

5%, 10-20%, and 30-50% centralities (2S1ET][^|^, and in the Pb-Pb collisions 
at y/SNN = 2.76 TeV with 10-20% and 40-50% centralities [TOl [U] [13]. The 
experimental data of pion momentum spectra, elliptic flow, and HBT radii impose 
very strict constraints on the parameters in the granular source model. They exhibit 
certain regularities for collision centrality and energy. 
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